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Abstract. The power spectrum of HI intensity fluctuation of the interstellar 
medium carries information of the turbulence dynamics therein. We present a 
method to estimate the power spectrum of HI intensity fluctuation using radio 
interferometric observations. The method involves correlating the visibilities in 
the u — v plane at different baselines. This method is particularly use full for eval- 
uating the power spectrum for the faint dwarf galaxies. We apply this method 
to 3 spiral galaxies and 5 dwarf galaxies. The measured power spectrum seem 
to follow a power law Pm{U) — AU", suggesting turbulence to be operational. 
Further, depending on the slope of the power spectrum, we expect the presence 
of 2D and 3D turbulence in those galaxies. 



1. Introduction 



Evidence has been mounting in recent years that turbulence plays an important 
role in the physics of the ISM as well as in governing star formation. It is 
believed that turbulence is responsible for generating t he hierarchy of structures 
present across a range of spatial scales in the ISM (e.g. 'Elmeg ren &: scalol 2004a : 
Elmegren scale 2004bl ) . In such models the ISM has a fractal structure and 
the power spectrum of intensity fluctuations is a power law, indicating that there 
is no preferred "cloud" size. 

On the observational front, power spectrum analysis of HI intensity fluctu- 
ations is an important technique to probe the structure of the neutral ISM in 
galaxies. The power spectra of the HI intensity fluctu ations in our own galaxy , 
the LMC and the SMC all show power law behaviour ( Crovisier &: Dickey 19831: 



GreenI Il993l : iDeshpande et al.l l2000l : lElmegreen et al . 
19991 ^ which is a characteristic of a turbulent medium. 



2001: Stanimirovic et al 



Recentlv lBegum et al. ( 20061 ) have presented a visibility based formalism for 
determining the power spectrum of HI intensity fluctuations in galaxies with ex- 
tremely weak emission. This fo rmalism was applie d to a dwarf galaxy, DDO 210 



and a spiral galaxy NGC 628 ( Dutta et al. 20081 ) . Interestingly, the HI power 



spectrum of both of these galaxies were found to be power law. In this report 
we present a elaborate description of this method. We also report the result of 
estimation of the power spectrum of 3 spiral galaxies and 5 dwarf galaxies. 
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2. A visibility based power spectrum estimator 

The specific intensity of the HI emission from a galaxy may be modelled as 

iu{e) = Wy0) [t + 5iv{0)\ ■ (1) 

Here we have assumed that this is the sum of a smooth component and a fluctu- 
ating component. We express the fluctuating component of the specific intensity 
as Wy{9) 5Iu{6), where 5Iy{6) is assumed to be a statistically homogeneous and 
isotropic stochastic fluctuation and W^{9) is the window function which quan- 
tify the overall large scale HI distribution of the galaxy. Since, the angular 
extent of the galaxies that we consider here is much smaller than the primary 
beam, we can write the visibility, the quantity directly measured by the radio 
interferometers, as 

V^(?7) = W{U% + W{U) 5X0) +Mu{U) (2) 

where the tilde ~ denotes the Fourier transform of the corresponding quantity and 
® denotes a convolution. In addition to the signal, each visibility also contains a 
system noise contribution Afu{U) which we have introduced in eq. ([2]). The noise 
in each visibility is a Gaussian random variable and the noise in the visibilities 
at two different baselines U and U is uncorrelated. At larger baselines, where 
the effect of the window function can be neglected, we can write, 

Vu{U) = WiU) SX{U) +M^{U) (3) 
We use the power spectrum of HI intensity fluctuations Pb_i{U) defined as 

{5X{U)6V{U')) = S\U - U') PniiU) (4) 

to quantify the statistical properties of the intensity fluctuations. We use angular 
averaging in place of the ensemble averaging denoted by the angular brackets. 
The square of the visibilities can, in principle, be used to estimate PmiU) 

{ V,{U)Vt{U) ) = \W,{U)\^ ® Phi{U) + (I M,{U) |2) (5) 

The last term (| J\fy{lJ) p), which is the noise variance, introduces a positive 
bias in estimating the power spectrum. The noise bias can be orders of magni- 
tude larger than the power spectrum for the faint external galaxies considered 
here. The problem of noise bias can be avoided by correlating visibilities at two 
different baselines for which the noise is expected to be uncorrelated. We define 
the power spectrum estimator 

= J (fU' Wy{U - U') W*{U + AU- U') Pm{U') (6) 

In our analysis we consider two different models for the window function of 
a galaxy. We present them here in the Fourier transformed form. 



(7) 
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and 

W,{U) = exp(-^2^2f^V2) (8) 

It is to note that the visibihties at two different baselines wih be correlated only 
if |Af7| < (vr^o)"^! and not beyond. In our analysis we restrict the difference 
in baselines to |A;7| < {^60)'^ so that W^{U + AU - U') ^ W^{U - U') and 
the estimator PHi(f^, AC/) no longer depends on AU. We then use the visibility 
correlation estimator 

PmiU) = {V,{U)Vl{U + AU)) 

= j cFU' I W,{U - U') p PuiiU') . (9) 

We use the real part of the estimator Phi(C/) to estimate power spectrum Piii{U). 
The real part is the power spectrum of HI intensity fluctuations convolved with 
the square of the window function. At large baselines, the effect of the window 
function can be neglected and we then have 

f>Hi{U) = C PuiiU) (10) 

where C = J \ Wiy{U) |^ (PU is a constant. The estimator PHi(f^) also has 
a small imaginary part that arises from the HI power spectrum because the 
assumption that Wu{U + AU -U) Wu{U — U) is not strictly vahd. We use 
the requirement that the imaginary part of Phi(C/) should be small compared 
to the real part as a self-consistency check to determine the range of validity of 
our formalism. 

The real and imaginary parts of the measured value of the estimator Phi(C/) 
both have uncertainties arising from (1.) the sample variance and (2.) the system 
noise. We add both these contrib utions to determ ine the 1 — a error-bars. The 
reader is referred to Section 3.2 of Ali et al.l ( 20081 ) for further details of the error 
estimation. 



3. Result and Discussion 

The dwarf galaxy data used here is from Giant Metrewave Radio Telescope 
(GMRT) observations. We use Very Large Array(VLA) archival data, for the 
spiral galaxies NGC 628, NGC 24 03 and NGC 1058. Details can be found in 
iDutta et al. 200^ ( in preperationi ) All the data are reduced in the usual way 



using standard tasks in classic AIPS Q- For each galaxy, after calibration the 
frequency channels with HI emission were identified and a continuum image was 
made by combining the line free channels. The continuum was hence subtracted 
from the data in the u — v plane using the AIPS task UVSUB. We correlate 
the visibilities in each channel as in eq. Q and then average over channels to 
get Phi(C/). We then fit a power law A U" to the estimator Phi(C/). The best 



^NRAO Astrophysical Image Processing System, a commonly used software for radio data pro- 
cessing. 
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Figure 1. Real part of the estimator, imaginary part of the estimator for 
the line channels and the real part of the estimator for the line-free channels 
are plotted from top to bottom respectively for the dwarf galaxy NGC 3741 
(left) and the spiral galaxy NGC 1058 (right). 1 — a error bars are also shown 
for the real part of the estimator for the line channels. Best fit is shown 
using a dark black line. Note that the best fit power spectrum for the galaxy 
NGC 1058 is a broken power law. 



fit A and a were determined through a minimization. To test whether the 
impact of the window function is actually small we have convolved the best fit 
power spectrum with | W,y{U) p. We estimate the goodness of fit (x^) to the 
data for the convolved power spectrum. The fit is accepted only after ensuring 
that the effect of the convolution can actually be ignored. The results are listed 
in table 1. 

The galaxies in our sample have dichotomy in the slope a r anging from 



-2.6 to —1.1. We have proposed a possible explanation for this in iDutta et al 



(120081 ') ■ This was based on the fact that DDO 210, where the power spectrum was 
measured across length-scales 100 — 500 pc, had a slope of —2.6 while NGC 628, 
a nearly face-on galaxy where the power spectrum was measured across length- 
scales 0.8 — 8 kpc had a slope of —1.6. We have interpreted the former as three 
dimensional (3D) turbulence operational at small scales whereas the latter was 
interpreted as two dimensional (2D) turbulence in the plane of the galactic 
disk. For a nearly face-on disk galaxy we expect the transition from 2D to 3D 
turbulence to be seen at a length-scale corresponding to the scale height of the 
galaxy. Continuing with this interpretation implies that we have also measured 
3D turbulence in NGC 3741, and 2D turbulence in NGC 2403, UGC 4459, GR 8 
and AND IV. 

The slope a = —1.0 ± 0.2 gives a good fit to the power spectrum measured 
at the length-scales 1.5 — 10.0 kpc for NGC 1058, where for the length scales 600 
pc— 1.5 kpc the best fit slope is a = —2.5 it 0.6. This, following the argument 
of dimensionality, is a transition from 2D to 3D turbulence. Since the length- 
scale 10.0 kpc is definitely larger than the scale height of any spiral galaxy, we 
interpret the wave-length —1.5 it 0.3 kpc at this transition as the scale-height 
of the HI disk of NGC 1058. To our knowledge this is the first observational 
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Dwarf 


DDO 210 


NGC 3741 


UGC 4459 


GR 8 


AND IV 


(1 a) a 

(2 a) range (kpc) 

(3 a) x'h^ 


-2.6 ±0.6 
0.10-0.50 

0.6 


-2.2 ±0.4 
0.15 - 3.75 

0.4 


-1.8 ±0.6 
0.16-1.78 

0.8 


-1.1 ±0.4 
0.1 - 1.5 

0.8 


-1.3 ±0.3 
0.56 - 6.2 

0. 1 


Spiral 


NGC 628 


NGC 2403 


NGC 1058 


NGC 1058 




(1 a) a 

(2 a) range (kpc) 
(3 a) x^lv 


-1.7 ±0.2 
1.0 - 10.0 
0.2 


-1.8 ±0.2 
0.5 - 4.4 
0.9 


-1.0 ±0.2 
1.5 - 10.0 
0.9 


-2.5 ±0.6 
0.6-1.5 
0.5 





Tabic 1. The results for the 5 dwarf and 3 spiral galaxies in our sample. 
Rows 1-3 gives (1) the best fit slope a, (2) length-scales over which the power 
law fit is valid and (3) the goodness of fit per degree of freedom. 



determination of the scale-height of a nearly face on spiral galaxy through its 
HI power spectrum. 

Wc have plan to use the same estimator to estimate the power spectrum of a 
large sample of spiral galaxy and to find the correlation of the slope of the power 
spectrum of those galaxies with the other measurable dynamical parameters like 
scale height, total HI mass etc. of the galaxies. 
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